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ABSTRACT 

The time delay between the formation of a population of stars and the onset of type la 
supernovae (SNe la) sets important Umits on the masses and nature of SN la progenitors. 
Here we use a new observational technique to measure this time delay by comparing the 
spatial distributions of SNe la to their local environments. Previous work attempted such 
analyses encompassing the entire host of cac;h SN la, yeilding inconclusive results. Our 
approach confines the analysis only to the relevant portions of the hosts, allowing us to 
show that even so-called "prompt" SNe la that trace star-formation on cosmic timescales 
exhibit a significant delay time of 200-500 million years. This implies that either the majority 
of la companion stars have main-sequence masses less than 3 Mq , or that most SNe la arise 
from double-white dwarf binaries. Our results are also consistent with a SNe la rate that 
traces the white dwarf formation rate, scaled by a fixed efiiciency factor. 
Subject headings: stars: evolution - supernovae: general 



1. INTRODUCTION 

Type la supernovae serve as the primary cos- 
mological standard candles (Colgate 1979; Branch 

& Tammann 1992; PhiUips 1993), due to their 
extremely regular light curves (Pskovskii 1977; 
Phillips 1993). However, this regularity need not 
imply that they share a common progenitor. In- 
deed, one striking bimodality is that the bright- 
est SNe la occur in star-forming galaxies, and the 
dimmest SNe la occur in galaxies with little star- 
formation (Hamuy et al. 1996; Howell 2001; van 
den Bergh & Filippenko 2005). This points to an 
evolution of SN la progenitors, and in fact, while 
several parameterizations of the SNe la rate exist 
(Greggio & Renzini 1983; Kobayashi et al. 2000; 
Greggio 2005), current data is well fit by a two- 
component parameterization (Mannucci et al. 2005; 
Scannapieco & Bildsten 2005), which takes the form 

SNRi^it) = AM4t) + BM,{t). (1) 

Here, we refer to the A-component, which is pro- 
portional to the total stellar mass of the host, as 
delayed., and to the B-component, which is is pro- 
portional to the instantaneous star-formation rate, 
as prompt. Note that these words have been oc- 
casionally used slightly differently in the literature 
(Greggio & Renzini 1983; Mannucci et al. 2006), as- 
signing SNe la as prompt or delayed/tardy based 
on timescales, rather than the parameterization in 
eq. (1). 

Despite the usefulness of eq. (1) on long time 
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scales, there are several reasons to believe that the 
prompt component should not be pinned to the in- 
stantaneous star-fomation rate. If core-collapse su- 
pernovae (SNcc) and prompt SNe la occured simul- 
taneously, the oldest stars in the Milky Way would 
be much more strongly enriched with iron group 
elements than observed (Mc Williams 1992; Scanna- 
pieco & Bildsten 2005). Furthermore, unlike SNcc, 
SNe la can not arise before stars evolve to form 
white-dwarfs, which takes at least 40 Myrs. Thus, 
it is likely that the prompt component exhibits its 
own delay time, r, such that 

SNRi^{t) = AM, (t) + BM, (t-r), (2) 

where current constraints, which use the global 
properties of host galaxies, place an indirect upper 
limit of T < 500 Myr on this timescale (Dahlen 
et aZ. 2004; Gal- Yam & Maoz 2004; Mannucci et 
al. 2005; Scannapieco & Bildsten 2005; Sullivan et 
oZ. 2006; Dilday et aZ. 2008). Galaxies older than 
this age have fewer SNe la, suggesting the charac- 
teristic delay time of the prompt component cannot 
be longer than 500 Myr. An alternative model is 
that SNe la occur at a rate that is directly propor- 
tional to the white dwarf formation rate (WDFR). 
In particular, Pritchet et al. (2008) argue for a uni- 
form conversion of newly formed white dwarfs into 
type SNe la at an efficiency of «1%. This would 
fix the la rate as being proportional to the stellar 
death rate of stars less than 8 Mq . 

In this paper, we employ a new observational 
technique for measuring the spatial distributions of 
SNe la and compare the results to analytical mod- 
els with the goal of constraining the delay time, r, 
for the B-component in the A-|-B formalism, and 
the feasability of the white dwarf formation rate 
(WDFR) model discussed below. The structure of 
this work is as follows. In §2, we review previous 
work constraining the progenitors of other classes 
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of transients, and build on these techniques to de- 
velop a new analysis applicable to type la SNe. In 
§3, we review the construction of the analytical host 
model for comparison from Raskin et al. (2008), in- 
cluding modifications necessary to make accurate 
comparisons to the results of our new observational 
technique. We present our results of the analysis 
and comparisons in §4, and conclusions are given in 
§5. 

2. DATA AND ANALYSIS 

2.1. Previous work 

Fruchter et al. (2006) developed an observational 
approach in which, for each host galaxy, they com- 
puted the fraction of the total light (or photon 
counts) in all pixels fainter than the pixel contain- 
ing a transient. By compiling a sample of such 
measurements into a cumulative histogram, they 
demonstrated that long-duration gamma-ray bursts 
are more likely than SNcc to cluster in the brighter 
regions of a galaxy. Kelly et al. (2008) expanded 
upon this analysis to distinguish between SNcc sub- 
types Ic, lb, and type II (see also Anderson & James 
2008). 

However, this study showed no difference between 
SNcc and SNe la, even though SNcc arise from mas- 
sive, short-lived stars rather than from white dwarfs 
that form from longer- lived stars. In fact, both pop- 
ulations followed the same distribution as that of 
the g-band light. In the core-collapse case, this is 
because young stars lead both to supernovae and 
the brightest regions. On the other hand, SNe la 
are not likely to arise from ongoing star formation, 
but they nevertheless share the same radial expo- 
nential profile as the g-hand light, which is caused 
by a radially-decreasing density of stars. This ra- 
dial gradient is large enough to obscure any signal 
caused by a delay in the prompt component (Raskin 
et al. 2008). 

2.2. Doughnut Method 

What is needed is a procedure that correlates 
SNe la with the properties of nearby regions, rather 
than with the host as a whole. In a spiral host, 
the ideal method for constraining SN la progenitors 
would be to measure the relative brightness of pix- 
els within annuli. In this case, as the density wave 
of star-formation moves around each annulus, SNe 
la would appear behind it at a characteristic sur- 
face brightness determined by the level to which a 
stellar population fades away in the g-band before 
SNe la appear. The g-band is ideal for this analysis 
as it fades away on the order of 100 Myrs, a similar 
time as that of white dwarf formation. However, ob- 
servations are never perfect, and observing a single 
annulus of a spiral host is subject to complications 
such as spurs, knots, and gaps, as well as the fact 
that stars rarely follow circular orbits. 

We call our solution to this problem the dough- 
nut method, and it builds directly on the method 
described in Fruchter et al. (2006). The idea is to 
expand an annulus radially by some small but ap- 
preciable radius, so as to encompass enough of the 




Fig. 1. — A sample image showing a spiral host, NGC 
3892, with a doughnut constraint, having a width of one-half 
the scale radius of the galaxy. The highlighted region has 
a small radial gradient, but still samples an area in which 
stars have an appreciable spread in ages reflected in their 
g-band surface brightnesses. The cross marks the location 
of type la supernova 1963J. 

host's morphological peculiarities to have a good 
representative sample, yet narrow enough to rep- 
resent local variations in the host light. Figure 1 
illustrates this concept. 

2.3. Sample Selection 

Our sample was drawn from the the Padova- 
Asiago Supernova Catalogue (Barbon et al. 1999), 
by selecting those events occurring within z < 0.07 
spiral hosts with Sloan Digital Sky Survey (SDSS) 
g-hand images. While this yielded 98 SNe la, we 
also removed those that occurred within 2" of a 
foreground star or within a galaxy inclined more 
than 60° or involved in a merger. This resulted in 
50 usable images, and as a control sample, we also 
selected all SNcc meeting these criteria, which re- 
sulted in 74 usable images. 

Note that before the advent of CCD astronomy, it 
was difficult to spot a supernova against the back- 
ground light when it occurred near the bright center 
of its host galaxy (see Shaw 1979). This induces a 
bias in most SN catalogs, the effect of which is to 
count fewer supernovae at small radii than would be 
expected from the stellar distribution in the hosts 
(e.g. Howell et aZ. 2000). Figure [2] shows the radial 
distribution of SNe la in our sample, which indicates 
a minor deficiency at small radii when our sample 
is compared to an exponential density distribution. 
Thankfully, our doughnut method overcomes this 
bias since it relies only on the local environment 
of the transient. Constraining our sample, then, to 
only recently discovered SNe has little to no effect on 
our results. On the other hand, when global prop- 
erties are considered, as in a traditional Fruchter et 
al. (2006) analysis, this bias can have drastic effects 
on the conclusions about SNe la progenitors. 

Given the typical resolution of SDSS images and 
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Fig. 2. — A cumulative histogram of the number of SNe 
la vs. radial coordinate showing a deficit at small radii 
when compared to an exponential density distribution. Note 
that while the density of stars in a spiral host follows an 
exponential radial profile, the projected area scales as the 
square of the radius, resulting in SN counts at the very 
center. 



the distances to these hosts, we settled on a dough- 
nut width of 0.5ro, or 0.25ro on either side of the 
transient, where tq is the scale radius. Host de- 
projection was carried out using the 25th mag B- 
band isophote as measured by the Padova-Asiago 
group, and to obtain the weighted average bright- 
ness within the uncertainties in the supernova loca- 
tion, we applied a gaussian convolution with a = 
1.2". Figure [s] shows the cumulative supernova 
distribution vs. the cumulative light distribution. 
Unlike analyses that use the full galaxy light, our 
approach is able to distinguish between the SNcc 
and SNe la distributions with 99.6% confidence, 
as quantified by a Kolmogorov-Smirnov (KS) test. 
SNe la are clearly biased to fainter pixels, while the 
SNcc distribution closely follows the g-band light 
distribution. 

In order to interpret this separation of the SNe 
la from the SNcc data, we construct an analytical 
model of a spiral galaxy following that of Raskin 
et al. (2008) with the necessary modifications to ac- 
count for the doughnut method. Such a model al- 
lows us to age date the SNe la by reproducing spa- 
tial distributions for a range of SN la delay times. 

3. MODEL DISTRIBUTION 

In such a model, we define the g-band surface 
brightness at a radius r and angle at a given time, 
t, as the convolution 



dt'Lg{t-t')±{t' ,r,e). (3) 



where Lg{t) is the g-band luminosity per unit so- 
lar mass of a population stars with a total age 
of t (Raskin et aZ. 2008), and S(t,r, 6*) is the star- 
formation rate surface density as a function of time 
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Fig. 3. — A histogram of the number of SNe la (crosses) 
located in pixels brighter than a given fraction of the 
g-band host galaxy light contained within the surrounding 
"doughnut," which shows that SNe la tend to occur in 
dimmer regions. A histogram of SNcc (circles) is also 
plotted, which closely follows the g-band light distribution. 
The two supernova histograms are different at the 99.6% 
confidence level. 



and position. For a galaxy with two spiral arms, 
this is well-approximated by 



S(t, r, 9) = So exp 



^ (4) 

where fip, the pattern speed, is given in eq, (|5|, and 
ro is again the scale radius, whence. 



i-^]-n.4r), (5) 



n (r) = l ^ ^ ^^'^'^ 

* 1 1'o/lkpc r < Ikpc, 

with Wo being the circular speed at Ikpc. 

In our simulations, we used a pattern speed that 
rotates 27r radians in 300 Myr, although this speed 
was varied by a factor of two and ruled out as hav- 
ing a significant effect on the results. Our model 
is described in further detail in Raskin et al. (2008) 
but these details are relatively unimportant as it 
is the characteristic brightness of the local environ- 
ment, rather than the angular position of SNe la 
that determine our results. 

In order to compute the SN la distribution, we 
replaced Lg{t) with SNRis.{t) as given by eq. p|, 
where the prompt component is modeled by a finSe- 
width gaussian distribution of width a — O.lr Myr, 
centered on an adjustable delay time, r. Larger 
values for a were also attempted, but these did not 
alter the best-fit r, and instead only reduced the 
goodness of fit, with the exception of very large a 
values, which yielded poor fits for all r. Further- 
more, it is crucial that the number of prompt and 
delayed la SNe in our sample be adjustable rather 
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than just adopting the average low-redshift num- 
bers. This is both becuase our hosts are selected to 
be star forming, and because of a likely Malmquist 
bias due to the fact that prompt SN la are signifi- 
cantly brighter than delayed SNe la. Thus we are 
left with two independent parameters: r, the delay 
time, and ^delayed the fraction of delayed SNe la in 
our sample. 

We can estimate this fraction in two ways. First, 
our SN la hosts have a mean color of — if = 3.53. 
Within this color bin SNe la occur at roughly three 
times the rate seen in elliptical galaxies (Mannucci 
et o/. 2005), which would suggest that ^delayed ~ 
0.3. Second, Howell et al. (2007) established a strong 
correlation between this fraction and the distribu- 
tion of SN la stretches at varying redshifts. They 
found that long stretches are indicative of prompt 
SNe la, while short stretches indicate delayed SNe 
la. Figure |4] shows a histogram of the stretches for 
a subset oi^ur sample that are available in the lit- 
erature (dark gray, Conley et al. 2008) and the dis- 
tribution for the same redshift bin found by How- 
ell et al. (2007, light gray). Again, because our se- 
lection criteria favor star-forming galaxies, a resul- 
tant two-component fit using the values for mean 
stretches (s) and dispersions (tr) of each component 
from Howell et al. (2007) demonstrates a mix more 
heavily favoring the long-stretch, prompt compo- 
nent than would normally be the case in a unbiased 
sample of low redshift hosts. Integrating these gaus- 
sians yields an estimate of -Fdoiaycd = 0.56. For our 
full analysis then, we consider i^doiaycd to be an ad- 
justable parameter, where our galaxy colors suggest 
i^doiayod ~ 0.3 and the results of our stretch analysis 
provide a natural upper limit of -Fdciaycd < 0.6. 

For every choice of -fdciaycd and r, two images were 
produced and analyzed using the doughnut method: 
one representing the host surface brightness and 
the other representing the SN la probability den- 
sity. The radial light decay across the doughnut was 
also measured for each observed host. For the most 
part, ln[L(rsw-0.25ro)]-ln[L(rsjv+0.25ro)] « 0.5, 
where rsN is the radial coordinate of the SN, but 
a few supernovafi were found in regions where this 
relation did not strictly hold, and this variation was 
used as a small modification to our model in a Monte 
Carlo fashion. However, the model cannot fully ac- 
count for the numbers of SNe la observed in pix- 
els fainter than the 20th percentile. In the model 
galaxy, with a doughnut centered on the annulus 
containing the SN, the dimmest regions will always 
lie at the outermost edge of the doughnut where 
SN probability is zero. In real galaxies, dim regions 
may lie anywhere, and a small number of the SNe 
la will be found in these dim regions. 

For the white dwarf formation rate model ( WDFR 
model), SNRig,{t) was derived from the Chabrier 
IMF (Chabrier 2003) convolved with the main- 
sequence turn-off ages from Bruzual & Chariot 
(2003) times 10% to approximate stellar lifetimes, 
for all stars below 8M0(see Raskin et a/. 2008). 
Stars whose turn-off ages are older than our simu- 
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Fig. 4. — A histogram of the distribution of stretches in 
our sample (dark gray) showing a different distribution as 
compared to that of Howell et al. (2007) for the same redshift 
bin (light gray). Using s and a values for each component 
from Howell et al. (2007), a two-gaussian fit for the prompt 
and delayed components with = 12.18 is overlaid. 

lated galaxy (wlO Gyr) do not contribute to the sim- 
ulated SN distribution. Pritchet et al. (2008) sug- 
gest a uniform conversion of ^1% of all white dwarfs 
to SNe la. However, the Fruchtcr et al. (2006) ap- 
proach, and by extension our doughnut method, dis- 
regards the total count of SNe, instead being con- 
cerned only with the relative distribution. 

4. RESULTS 

In Figure [5] we plot A+B model distributions with 
-Fdciaycd = 0.3 and varying delay times as well as 
the result of the WDFR model. It is clear from 
the A-l-B model curves that larger values for r pro- 
vide increasingly good fits to the data, while the 
50 Myr model is inconsistent with the observations. 
Because our model galaxy has two free parameters, 
however, we must account for degeneracy between t 
and -Fdciaycd- Figure [6] illustrates the two-parameter 
KS probability map for our sample when compared 
to our model with varying -Fdciaycd values and delay 
times. 

Both the two-parameter map and the average KS 
probabilities for all tested -Fdciaycd values show a 
clear trend toward longer delay times with the prob- 
ability reaching its maximum at ~ 500 Myr and one 
value of Fdciaycd yielding a probability of 82% at 
this T. This is a significant delay time, roughly cor- 
responding to the lifetime of a 2.5 M© star at solar 
metallicity, and is much longer than the minimum 
theoretically expected time of 40 Myrs. The best fit 
delay time for all models, even those with high val- 
ues of i^deiayed, IS longer than 200 Myr and probably 
longer than 300 Myr. On the other hand, global 
approaches place an upper limit of r < 500 Myrs 
(Dahlen et al. 2004; Gal- Yam & Maoz 2004; Man- 
nucci et al. 2005; Scannapieco & Bildsten 2005; Sul- 
livan et al. 2006; Dilday et al. 2008). Relating a 200- 
500 Myr delay to main sequence lifetimes, it is clear 
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Fig. 5. — Model curves representing the A+B model for 
several delay times as well as the white dwarf formation rate 
model are shown compared to the observed distribution of 
SNe la. Both the r = 300 and WDFR models provide excel- 
lent fits to the data. 

that majority of SNe la stars have main-sequence 
masses less than 3 Mq , or alternatively, that most 
SNe la arise from double- white dwarf binaries. It 
is important to note that while our model provides 
strong constraints on the characteristic SNe la delay 
time, it does not rule out that a fraction of SNe la 
occur at shorter delays (Aubourg et al. 2008; Man- 
nucci et al. 2005; Anderson et al. 2007; Anderson & 
James 2008; Totani et a/. 2008). 

In Figure [5] we also plot the WDFR, whose sim- 
ilarity to the A-t-B model with t = 300 Myrs is 
striking. One cannot truly be ruled out in favor of 
the other. A KS test for the WDFR yields a value 
of 72%, and in fact, the average delay for this model 
as calculated by ]^ SNR{t)tdt/ SNR{t)dt is 
« 500 Myr. 

Finally, to place an upper limit on the frac- 
tion of prompt SNe la with very short delay times 
we considered a three-component model in which 
SNR{t) = AM.^{t)+BilvU{t-imMyvs)+B2M*{t- 
300Myrs). By fixing r, we again consider two free 
parameters, ^delayed and Bi/BTot, the fraction of 
the prompt component with a delay time of 100 
Myrs. For large values of Bi/Btoi, a KS test yields 
a low statistical likelihood in the model for all val- 
ues of -Fdciayod- Wc Can Say with some confidence 
that the maximum allowed value of Bi/Btoi is 0.3, 
corresponding to 30% of the prompt component to 
coming from a 100 Myr old population, while the 
parameter space with the highest probability has 
Bi/BTot = 0. 

5. CONCLUSION 

Differentiating the spatial distribution of SNe la 
from core-collapse SNe is a difficult problem, both 
due to the stochastic nature of SN la detections 
and the potential for the delayed component of the 
la rate to obscure the prompt SN la spatial dis- 
tribution. However, by restricting the Fruchter et 
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Fig. 6. — Top Panel: A KS probability map for the A-l-B 
model with varying parameters Fdolayod ^'nd t. Degeneracies 
between the two fit parameters are extremely small. For all 
values of Fdoiaycdi T 5: 100 Myrs is very unlikely, and the 
best overall fit value is 500 Myrs. The solid contours mark 
various levels of KS confidence. Bottom Panel: Average 
and best fit KS probabilities over a range of i^dolayod values, 
each as a function of r, again showing that substantial delay 
times are highly favored. 

al. (2006) analysis to an annulus, our "doughnut 
method" seperates the local host properties from 
global properties, allowing us to measure the im- 
pact of the delay time of the prompt component on 
the spatial distribution of SNe la. By comparing 
our observations to analytical models, we have es- 
tablished a strong case for a modified A+B model 
in which the prompt component is delayed by 300- 
500 Myrs. Note that this is an average time and 
that there can be a considerable spread in this value. 
Thus, a three-component model allows up to 30% of 
the prompt component to have shorter delay times 
« 100 Myrs, as suggested by Auborg et al. (2008) 
and Mannucci et al. (2006). Alternatively, we also 
found that a model in which the la rate is directly 
proportional to the white-dwarf formation rate, as 
considered by Pritchet, Howell & Sullivan (2008), 
reproduces the observed spatial distributions of SNe 
la very well. 

Regardless of which model proves to be the best 
fit, it is this characteristic timescale that is most 
important for calculations of cosmic enrichment, 
and the results of our analysis using both the A-l-B 
model and the WDFR model are evidence of a long 
characteristic timescale. One caveat is that all our 
measurements have been made in relatively high- 
met allicity grand-design sprial galaxies, and several 
observational and theoretical studies have hinted 
that SN la rates and properties may be substantially 
different at lower metallicities (Timmes et al. 2003; 
Gallagher et aZ. 2008; Cooper et aZ. 2009). There is 
still much more to be learned about these ubiqui- 
tous but mysterious cosmological probes. 
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